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Abstract 

We write a multiple integral formula for the partition function of the Z-invariant six vertex 
model and demonstrate how it can be specialised to compute the norm of Bethe vectors. We 
also discuss the possibility of computing three-point functions in A/" = 4 SYM using these 
integral formulas. 



^ ! 1 Introduction 

\ The six vertex model is one of the most studied solvable lattice model in statistical mechanics [1] . 

Using the Bethe ansatz for eigenvectors and eigenvalues of the transfer matrix, Lieb and 
Sutherland [1] solved the regular square lattice zero-field six-vertex model with periodic boundary 
conditions. By computing the largest eigenvalue they obtained explicit solutions for the free-energy 
per site in the thermodynamic limit. 

The Bethe ansatz solution for a periodic square lattice with sites requires that for each 
eigenvalue one has to solve n non-linear equations in n unknowns (the 'wave numbers' fci, . . . , fc„). 
' This can only be done explicitly in the large-L limit under the 'string hypothesis'. Usually one 

^ \ can solve those equations explicitly only for the largest and near-largest eigenvalues, which limits 

one to considering the full infinite square lattice. 

On the infinite square lattice vertex operator methods are available (see e.g. [5]), which have 
resulted in multiple integral expression for correlation functions of the model. 

The situation is different when one considers fixed boundary conditions. Izergin and Korepin 
constructed an explicit determinant solution for the case of 'domain wall boundary conditions' [5117] . 
Alternative representations for this same partition function have also been proposed by other 
authors. See for instance [8HT5] and references therein. 

In |14| . Baxter constructs the partition function for a general Z-invariant six- vertex model, 
which is given explicitly (apart from a simple factor) by a Bethe ansatz type expression. However, 
in this case the 'wave numbers' do not need to be evaluated from a complicated set of simultaneous 
equations: they are known explicitly. The disadvantage of Baxter's expression however is that it 
is given in terms of a large summation over the symmetric group. 

The main results of this paper is to rewrite Baxter's expression as a multiple contour integral 
over a factorised polynomial kernel. As an application of this result we specialise it to the regular 
square lattice with two specific boundary conditions. This allows us to write a new representation 
for the partition function of the six vertex model with domain wall boundary conditions as a 
multiple contour integral. Likewise we derive an off-shell multiple integral expression for the 
scalar product of two Bethe vectors of the six-vertex model transfer matrix. 



1 



Such formulaeare important as recent developments in the computation of three-point functions 
in A/" = 4 Super Yang-Mills (SYM) seem to require manageable expressions for the norms of Bethe 
vectors for general values of the parameters. With this in mind, in Section[B]we illustrate how the 
integral formulspresented here can be embedded in the framework of |15) . 



2 The six vertex model 



The exact solvability of a vertex model in the sense of Baxter 
solutions of the Yang-Baxter equation. This equation reads 



is intimately related to the 



(1) 



where Cij G End(Vi Cg) Vj). The complex variables A and fi correspond to spectral parameters 
while Vi denotes a complex vector space. We refer to the solutions of (H]) as /^-matrices and for 
the six vertex model the corresponding £-matrix is invariant under the C/5[s[(2)] algebra in the 
fundamental representation. In that case Yi ~ and the associated /^-matrix explicitly reads 



C = 



fa 



V 





a J 



where 



a(A) — sinh (A -I- 7), 6(A) = sinh (A) and c(A) = sinh (7). 



(2) 



(3) 



To characterise the vertex model associated with a solution of the Yang-Baxter equation we 
shall employ the following notation, 



^= ^ ^ij eafs <8) e^s 



(4) 



where Bap = \a) {I3\ where the vectors \a), as well as their duals, form a basis of V;. In this way 
we associate the Boltzmann weight £^1^ to the vertex configuration {a, f3,j,S} as shown in the 
Fig.ffl 



Figure 1: Vertex {a,(3,j,6} and its Boltzmann weight. 



For the six vertex model the indices a and /3 in Q run through the set {f, ],} while 7 and S 
run through {—>,•<—} defined as 

It) - K) = (J) e V, and 1;) = 1^) = (^^ e V. . (5) 

From ([2]) we thus have six possible vertex configurations which are depicted in the Fig. [2l The 
vertices described in Fig. [2] can be juxtaposed in a two-dimensional lattice with N rows and L 
columns as shown in the Fig. |31 The probability of having the vertex {aij,aij+i,f3ij,f3i+ij} at 
the Cartesian coordinates is weighted by the factor and the partition function 
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Figure 2: The 6 vertex configurations of the model. 

/5jv+i,i ... ... /3jv+i,i 



<^1,L + 1 



Figure 3: Two dimensional square lattice with boundary conditions index by aij and /J^.j who 
take values in {^,— >•} and {t,!} respectively. 



of the system is formed by the summation over all possible configurations of the product of vertex 
weights. More precisely, the partition function Z is defined by 



N L 
{a,,,,ft.,}i=lj = l 



(6) 



It is possible to generalise this construction in the following way to include inhomogeneities 
while maintaining integrability. To each horizontal line i of the square lattice we associate a 
variable A^, and to each vertical j a variable ^j. The local Boltzmann weight for a vertex 

{ttij-, ttij+i, is then given by £"''^.^3^^'^^'^ (Ai — /ij) and the inhomogeneous partition 
function is defined as 



N L 



(7) 



Boundary conditions play a fundamental role in the evaluation of the partition function 0. In 
this paper we will consider periodic boundary conditions in the horizontal and vertical directions 
(c^i.L+i = CKi.i and /Jat+i j — as well as fixed boundary conditions with specified values of 
aij and f3ij at the borders. 

2.1 Monodromy and transfer matrix. 

Each row of the two-dimensional lattice formed by the juxtaposition of vertices can be conveniently 
characterised by a matrix 7~(A) = T(A, {/ij}) usually referred to as monodromy matrix. This 
matrix has components 
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Figure 4: Diagrammatic representation of the monodromy matrix. 

which is diagrammatically represented in the Fig. 21 In this way the two-dimensional lattice 
described in the Fig. [3] can be completely built in terms of the operators T(A, 
In a more compact form we have 



^(A, {^ij}) = £^i(A - ^i) . . . Cal{\ - Mi)- 



(9) 



which is an operator on the tensor product space (g) Vi (g) • • • (8) Vl. The space is usually 
called the auxiliary space while we shall refer to the tensor product Vi (g • • • (SD Vi, as the quantum 
space. Thus the monodromy matrix T{X) is a matrix in the auxiliary space whose elements 7^* 
act on the quantum space. In the diagrammatic representation of the monodromy matrix given 
in Fig. m the contraction of indices labelling the horizontal edges represent the product in the 
auxiliary space while the product in the quantum space is given by the contraction of indices 
labelling the vertical edges. For the six vertex model the monodromy matrix can be conveniently 
denoted as 

(10) 



Here and in the following we suppress where possible the dependence on {/^j} for clarity of notation. 



3 Periodic boundary conditions 

The evaluation of the partition function ([6]) depends drastically on the boundary conditions chosen. 
In this section we shall consider periodic, or toroidal boundary conditions (a^^L+i = a^.i and 
Pn+i.j — and review the basics of the quantum inverse scattering method (QISM) for 

computing the six- vertex partition function in this case. 

From ^ and it is a standard construction that the partition function for toroidal boundary 
conditions can be written as 



{^J■J}) = Trc2®t (r(Ai) • • • t{Xn)) , 
where t(A) is the transfer matrix defined by 

t(A) = t(A, {a*,}) = Tr^r(A) - A{X) + D{X). 



(11) 



(12) 



Due to the Yang-Baxter equation the transfer matrices T(Ai) all commute and are simultaneously 
diagonalisable. The partition function is therefore given by 



2^ N 



z{{x,},{^^J})^J2]lM^^), 



(13) 



where is the fth eigenvalue of the transfer matrix r. 
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3.1 Bethe ansatz diagonalisation 

Let us first define tfie 'pseudo-vacuum' states and jifj^) by 

L L 

m = (g) It) and \^{) - (g) li) . (14) 

i=l 1=1 

An important ingredient in tfie QISM setup for diagonalising tfie transfer matrix r(A) is tlie action 
of tlie monodromy matrix elements on tlie states Y^^) and |!f'4,). For \^^) tliese are given by 

L 



A{X) \^^) = \[a{\^ f,,) B{X) \^^) - * 

L 

C(A) Itf't) = D{X) Itf't) = n - ^^■) I'^t) • (15) 

The relations follow from the definitions ^ and PH)) together with the triangularity of C^j 
on local states. 

Due to the Yang-Baxter equation ([1]) one can readily show that the monodromy matrix ^ 
satisfies the following quadratic relation 

7^(A - jy) [r(A) Tii^)] = [r{v) ® r(A)] nix - v) (le) 

which is known commonly as the Yang-Baxter algebra. Here TZ = PC where P stands for the 
standard permutation matrix. Thus (jl6p yields commutation relations for the elements of the 
monodromy matrix (1101) . Among the relations encoded in (|16p we shall make use of the following 
ones, 

A{X)B{.) = "^^^B{.)A{X) - ±-^^B{X)A{.l 
B{X)B(y) = B(y)B{X) . 

We are now in a position to state the main theorem of the algebraic Bethe ansatz, or quantum 
inverse scattering method. 

Theorem 3.1. Let the numbers Ai, . . . , A„ be solutions of the equations 
Then, 

\Xi,...,X„) = B{Xi)---B{K,)\W^) (19) 
are eigenvectors of the transfer matrix r(A, {/ij}) with eigenvalues A{X, {fij}) given by 

MK M) = n «(A - n 't!^ + n ha - „, n . im 

3 = 1 k=l ^ ' 3 = 1 k=l ^ ' 



Proof. The theorem follows from the definition (112]) of the transfer matrix, the action dTSl) on the 
pseudo- vacuum, and relations ()17p . In particular, the expression for the eigenvalue originates from 
([T5|) and the first terms on the right hand side of (|17p . Equations (fTSl) imply that unwanted terms 
arising from the other terms on the right hand side of ([TT]) cancel. □ 
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An important problem in the theory of solvable lattice models is to provide manageable ex- 
pressions for correlation functions. These are of the form 

(j.i,...,z/„|0|Ai,...,A„), (21) 

where O is some operator. The simplest case concerns an expression for the norm of Bethe states, 

^({i^.},{A,}) = (j.i,...,z.„|Ai,...,A„). (22) 

Below we shall see that ((22|) is equal to a special case of the partition function ([7]) with fixed 
boundary conditions. In the next section we show that such partition functions can be realised as 
multiple contour integrals over a factorised kernel. 



4 Fixed boundary conditions 

In this section we review some known results for the partition function ^ in the case of fixed 
boundary conditions. In particular we recall Baxter's formula for such partition functions for 
general graphs and boundary conditions, expressed in terms of a sum over the symmetric group 

M- 

Let us first look at the well known case of domain wall boundary conditions. 



4.1 Domain wall boundary conditions 

For the case of domain wall boundaries with N = L, the edges on the vertical and horizontal 
borders in Figure [3] assume a particular configuration. Naturally one needs to be careful when 
choosing the boundaries since an inappropriate choice can render a trivial (null) partition function 
or force the system to freeze in a particular configuration. Domain wall boundaries for the six 
vertex model were first introduced by Korepin [7j, and in our notation correspond to the choice 
ai^i =— i>, I^L+ij =t and ai^L+i =■<— , =4-- This special boundary condition is illustrated in 
Fig. [SI 



Figure 5: Domain wall boundary condition. 

Using the prescriptions given in Sec.[2l we can see that the top and bottom boundary represent 
the states |tf't) and defined in (fT4|). The left and right boundaries in Figure [5] select the ele- 
ments of the monodromy matrix forming the operator O = Jli^i B{Xi, {Mj})- Thus our partition 
function can be written as 

L 

^({A.}, {a*.}) = (f'il n ^(^" iA*.}) I'^t) , (23) 

i=l 
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where we emphasise that the rapidities {A^} may be chosen arbitrarily and do not have to satisfy 
the Bethe ansatz equations. It is weh known that (j23|) can be written as the Izergin-Korepin 
determinant [BJ[7] as well as sums over the permutation group [TUlll3j . 

4.2 Scalar product 

The scalar product (1221) can be written in a similar fashion as the domain boundary partition 
function. From the definition of the Bethe eigenstates in (fT9|) it follows immediately that the 
scalar product (|22|) can be written as 

n n 

N{w^}, {A,}) = (fti n ^("^^ {^fc}) n ^(^^■' i-^t) > (24) 

i=l 1=1 

where the sets of numbers {vi} and each satisfy the Bethe equations p^ . Graphically (|24l) 
is depicted in Fig. [51 

1^1 ... f^L 



■< K 



-< Ai 



A, A. 



Figure 6: Scalar product boundary condition. 



4.3 Perimeter Bethe ansatz 

In this section, we review the result of Baxter [l^ for a general Z-invariant six-vertex model, 
adjusting some of his notation to serve our purposes 0. The general model is considered on a 
simply connected convex planar domain D, with m straight lines within it, starting and ending at 
the boundary [H1IT7] . No three lines are allowed to intersect at a common point, (see Figure [7] for 
an example). We fix some point on the boundary, and label the ends of the lines as 1,2,..., 2m, 
in a clockwise direction. The line with endpoints i and j is referred to as 'the line where 
1 !i * < j ^ 2m. We associate a 'rapidity' Vi with each endpoint i, and for each line («, j) impose 
the constraint 

Vj = Vi - 7. (25) 

^in particular, note that Baxter labels the endpoints in an anti-clockwise direction, and writes his formula with 
respect to the left rapidities 
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Figure 7: Example of a general planar graph inside a domain T> without multiple crossings at one 
point for to = 6. 

Baxter uses the terminology 'right' rapidity and 'left' rapidity of a line («, j) for Vi and Vj respec- 
tively. We also define the set of all rapidities V — {ui, . . . ,W2m}, as well as the set U of all right 
rapidities. 

A six-vertex model is constructed by placing arrows on the edges of the lattice so that each 
vertex (intersection) has two in-pointing arrows in and two out-pointing arrows (the ice-rule). A 
vertex at the crossing of lines («, j) and (fc, I) with i < fc < j < as in Figure [8l is given a weight 
w{i,j\k,l) defined below. The six possible arrangements as viewed from the boundary points k 



I 




k 



Figure 8: A vertex at the crossing of lines («, j) and (fc, /) with i < k < j < I has weight w{i,j\k, I). 
and j are shown in Figure [9] Boltzmann weights are assigned to the six configurations as follows 

P' J> 

n n n n ^ ^ 

Wi W2 W3 W4 W5 We 

Figure 9: Six vertex configurations on vertices inside V. 
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(in the notation of this paper): 

Wi = W2 ^ 1, 



H^i - Vk) sinh(w; - Vk) 

Wi=W4. = — = — ■ — (26) 

a{vi - Vk) sinh(u/ - Vk + 7) 

c{vi - Vk) sinh(7) 
wz=Wg- - 



a{vi - Vk) sinh(w/ - Wfe + 7) 

We also fix boundary conditions, with m arrows pointing into I) and m arrows pointing out. 
We specify these boundary arrows with the set of endpoint locations where an out arrow occurs, 
X = {xi, . . . ,Xm}- The partition function of the model is then a function of X and V. We now 
come to the main result (for our purposes) of 14 : 

Theorem 4.1 (Baxter). Let V — {vi, . . . ,W2m} be the set of all rapidities, and let U be the set of 
right rapidities. The partition function Z{X\V) for the Z invariant six-vertex model is given by 

ZiX\V):^J2 n ^(*'J>'0 

cfgs. {i,j\k,l) 

/ \ m 

=^-^E n i^^n^^-(-o (27) 

P l<i<j<m ^ ^ ^' i=l 

where 

x — l 2m 

<f>xiu) = a{vj -u) Y\_ b{vj - u), (28) 

and the pre-factor is defined by 

C= (-l)"tsinh(7)-™ J] a{u,-uj)^, (29) 

l<i,j<m 

where n-^ is the number of out-point arrows with coordinate x < m and the variables Ui refer to 
the elements of the set of right rapidities U . 



5 Integral expression 

In this section we will derive a multiple integral expression for the partition function ([7|) in the case 
of fixed boundary conditions. Theorem 14.11 provides an explicit expression derived by Baxter |14) 
for such partition functions for general graphs and boundary conditions. This expression is in 
terms of a large sum over the symmetric group and in this form not immediately useful. We hope 
our integral expressions are more manageable and may lead to further progress for correlation 
functions as well as higher rank cases. 

We wish to write the sum over permutations as an integral, summing over residues. We 
introduce auxiliary integration variables Wi, ? = 1, . . . , m. This leads to the following result. 

Theorem 5.1. Let V = {fi,...,W2m} be the set of all rapidities, and let U be the set of left 
rapidities. 

7fYim / / n" 1 n,"L.+i «(^.- ^ ^^)b{m w,) ^ ^ 

= J2^f---f n^^'^"^^ d^i...du;„. (30) 

The contours of integration are around the poles at Wi — uj, i,j = 1, . . . ,m and the prefactor is 
again given by 

C= (-l)"^sinh(7)-" H a{u,-uj)\ (31) 

l<i,j<m 

and 4>xiw) and are as in Theorem \^.l\ 
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Proof. The factor Y[j>i ^{'^i i'^ the denominator imphes that the only nonzero contributions 

arise from residues at different poles, i.e. if {wi, . . . , Wm} is a permutation of {ui, . . . , Um\- □ 

5.1 Domain wall boundary conditions 

We now consider the result of Baxter il4', specialised to a L by L square lattice with domain wall 
boundary conditions. Using the notation pLj = we have AL rapidities though only 2L are 

independent. The full set of rapidities are given by 

V = {^1, . . . Ai, . . . , AL,/ii - 7, ■ • ■ - 7, Ai - 7, . . ., Al - 7} (32) 

while the set of left rapidities is 

U = {wi, . . . , u2l} = {/xi, . . . ,^L, Ai, . . . , Al}. (33) 

We take boundary conditions of out arrows being on the vertical lines, and in arrows on the 
horizontal lines according to Fig. [H Thus the positions of the out arrows are 

X = {l,...,i,2i+ l,...,3i}. (34) 

With these definitions we have that the domain wall boundary partition function ()23|) is equal to 



L 

z{{K}, {^^,]) - n - ^^j)z{x\v). (35) 



The domain wall boundaries specialisation of (j30l) implies that the only nonzero contributions 
arise from residues at Wi = /i^ for i = 1, . . . , L, and if • • ■ , W2l} is a permutation of 

{Ai, . . . , Al}. With this specialisation, and the re-labelling of the integration variables WL+i Wi, 
the formula becomes 



^({AJ,{/i,}) 



/ sinh7\ ^ / / ULi Uj=^+l aiwj - Wi)b{wj - Wi) 



L 

3 = 
L i-1 



^W_W_b{^j ~ Wi) J]^ a{'Wi - ^ij) dwi . . .Awl (36) 

i=l 3 = 1 3=1+1 



where the integrals are around the poles at = A 



5.1.1 Homogeneous limit 

The integral formula enables us to readily calculate the homogeneous limit of the partition 
function. With the changes of variables 

1, f\-tx\ , 1, fl-tz\ 1, fl-ty\ 

--21°HT^)' '=21°HT^)' ^=2l°HT^)' ^''^ 
we obtain 



Z 



tL 

L L 



iVj - Zi)(l + KZi + yjZi) 



X 
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where n — — (t + 

We now set Zi — z and yi = y for all i (i.e. Xi = X, fj,i = fi), 



T ^ 



L 



(?/ - z)^ (1 + + yz)^ -j-j-^ 



(1 + Ky + y2)L/2(l + KZ + z2)-£'/2 J^i (Zj- - Xi) 

, i A 1 

(1 + nxi + x,y)^ ^J: [y - x,)^+ 



X TT -; h-TTt -i—r-dxi...dxL- (39) 

i-^ l + KXi+x,y)^ i-i-. y-x,)L+^-' ^ ' 



Also the limit a = h = c = 1, where the partition function counts the number of alternating sign 
matrices [TH], is obtained as follows. We first set a(A, ^) = h{X, i.e. X — ^ = —7/2 + i7r/2. In 
the transformed variables we achieve this by setting z = and y — 1, corresponding to A = —7/2 
and /i — — i7r/2. We also normalise each of these weights to 1 by dividing through by a(A,/i) and 
thus obtain 

L 



(2 + «)^V2(2,ri)^7---7 nti(-^.)' 

L 

[](1 + (k + - dxi . . . dxL. (40) 



X 

1=1 



(Note that we would finally get (a = 6 = c = 1) by setting 7 = Itt/S, i.e. k = —1). Alternatively 
we set K = r — 2, and write equation PO)) as the following constant term expression: 



Ai(r) = CT 



n n (a^i -a;z)(l + (r-2)xj +Xja;j) 



n 



i=l -^J 



{1 - x,)^--+^i + {t - i)x,y 



, (41) 



which yields polynomials in r corresponding to the enumeration of alternating sign matrices, 
by number of —Is. This latter formula is very similar to expressions obtained in |19H23| as 
polynomial solutions of the g-deformed Knizhnik-Zamolodchikov equation and their relation to 
the combinatorics of alternating sign matrices and symmetric plane partitions. 

5.2 Scalar Product 

Next we consider the result of Baxter [14] specialised to the 2nx L rectangular lattice depicted in 
Fig- El This specialisation renders the norm of Bethe vectors (|24|) . and according to the conventions 
of Section we set m = i + 2n. We have rapidities 



V = {/ii, . . . Ai, . . . , A„,i^i, . . . ,z/„,/^i - 7, . . . - 7, i/i - 7, . . . ,i/„ - 7, Ai - 7, . . . , A„ - 7} 

(42) 

and the set of left rapidities is then given by 

U = {^i,...,/iL,Ai,...,A„,z/i,...,t'„} . (43) 

Here we use the definitions ftj = ^L+i-j, Xj = Xn+i-j a-nd Vj = Vn+i-j. In order to characterise 
the scalar product p4| we also need to specify the location of the out arrows X. We have these 
out arrows appearing in three distinct blocks with positions: 

a: {1, . . . , L, L + n + 1, . . . , i + 2n, 2L + 2n + 1, . . . , 2L + 3n} (44) 
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Similarly to the case of domain wall boundaries, the norm (j24p consists of the partition function 
Z{X\V) up to a normalisation factor arising from the normalisation of the statistical weights. Thus 
we have 

n L 

N{{l^^}, {A,}) = n n «(^'' - ^^X^' - A^^)^(^l^)- (45) 

with Z{X\V) computed using (gl]), (gH]) and (gH). With this specialisation the integrand of ((50)) 
has poles only at Wi = fjn for i = (see Appendix VK\ and thus the integration over 

this particular subset of variables is trivial. After performing the trivial integrations and some 
re-arrangments detailed in Appendix]^ we are left with 



.r/r . rx ^^ (-l)(^-')"(sinh7)2" " 1 



(27ri)2" .11^ a{\i - Vj)a{vj - \i)a{v, - v^Y 

n?=l Y{%^+l a{Wj - Wi)h{Wj - Wi) 



dwi . . . Aw 



2n- 



n L 71 n 

^WW a{m - Hj)b{^jij - Wn+i) W W a{wi - i'j)a{yj - Wn+i) 

n n 

xY\_ Y\. ^^^^ ~ Wi)b{w2n+l-i - Vj)a{vi - Wj)a{w2n+l-j - Vi) . (46) 
i—1 



6 Three-point functions in = 4 Super Yang-Mills 

The formulas (gni) for the scalar product of Bethe vectors is valid for arbitrary values of the com- 
plex parameters {I'i} and {Xj}. On the other hand, the expression (|46p will yield the norm of the 
transfer matrix eigenstate (ITOl) only when the set {vi} is related to the set {Xj}, i.e. i/j = A*, and 
with the set {Xj} subjected to the Bethe ansatz equations ([T8)). In that case Slavnov's formula [24] 
provides a rather simple expression for the scalar product in terms of a determinant of a n x n 
matrix. For generic parameters {vi} and {A^}, one would have to consider the formulsof [71125] . 
given in terms of a large sum over partitions. However, the recent developments in the compu- 
tation of three-point functions in A/" = 4 Super Yang- Mills (SYM) seems to require manageable 
expressions for the norms of Bethe vectors for general values of the parameters and here we shall 
illustrate how the integral formula (|46p can be embedded in the framework of [15] . 

Three-point correlation functions of primary operators in A/" = 4 SYM theory are constrained 
by conformal invariance to have the form 



where J\fj are normalization factors and Ai are the respective conformal dimensions. The di- 
mensions Ai are already fixed in the computation of two-point functions and a large literature is 
devoted to that problem in the planar limit of the A/" = 4 SYM, exploring integrable structures. 
For a detailed revision on this subject see for instance |2B] and references therein. Here the quan- 
tity we wish to evaluate is the structure constant Cijk, which has been worked out to leading order 
in [15] employing the spin chain picture of single trace operators introduced in [27] . See also [28] 
and [29] for previous works on A/" = 4 SYM three-point functions within the spin chain picture. 

The structure constant Cijk in the Af = A SYM admits the perturbative expansion ajk — 
W'^ijk + W'^^iji ~^ W'^^iji + ■ • where and 6 are respectively the number of colours and the 't 
Hooft coupling constant. We shall consider the leading order term c^js described in [15] and in 
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order to compute it from our results we first need to consider a special limit usually refereed as 
rational limit. We start by performing the change of variables 



7 = 17 , 



A — 7(it — ^) and 

and we also set the inhomogeneities fij =0. 
From ©-(Uni) we find that 

lim . }_. Ti\) - Tiu) , 



+0 sinh7' 
-+o ' 



(48) 



(49) 



where T{u) corresponds to the monodromy matrix considered in |15| . In practice the above limit 
corresponds to considering our results with statistical weights a, b and c given by 



a = 
b = 
c = 



(50) 



Now we consider three independent XXX-type spin chains with periodic boundary conditions 
which we shall refer as Ai, A2 and ^3. Those spin chains can be obtained as the logarithmic 
derivative of the transfer matrix (jl2p up to normalisation factors after considering the rational 
limit above described. This is a standard construction and a detailed description can be found 
in [2F. Within the spin chain picture introduced in [27], each single trace operator Oi in the 
SU{2) sector of the A/" = 4 SYM at one-loop order is then represented by the spin chain Ai under 
the zero momentum condition. 

The state of each spin chain Aj with Lj sites is characterised by rapidities satisfying Bethe 
ansatz equations. For instance we use the set of mi rapidities {ui, . . . ,Umi} to characterise the 



spin chain Ai, the m2 rapidities {vi, 



,} for A2 and 7713 rapidities {wi, 



J for A3. 



They are subjected to the following Bethe ansatz equations: 



b{ui) 



Li 



1 L2 



b{wi) 



n 

fc=i 

1T12 

n 
n 

k^i 



a{ui 



Uk 



a{uk 


- Ui -\ 




a{vi 


-Vk^ 




a{vk 






a{wi 


- Wk 





a{wk -Wi + ^) 



i = 1, . . . ,m\ 



. ,m3 



(51) 



and the zero momentum condition translates into 



mi , s 



™2 / N 



n 



a{wi) 



(52) 



We remark that now we are considering the statistical weights a, b and c given in (|50p. Next we 
define the functions 



E2^\{ 



a{vj) 



3 = i 

m-i mi 



a[ui - Uj + ^) 



^i=n n 



F2 



G2 



"12 "12 / 1 i ^ 

TT TT - ^3 + 2) 

"i2 "12 / I i ^ 

-TT -TT - + 2) 



(53) 
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and separate the set of rapidities {u} — {ui, . . . ,Umi} into two subsets {a} and {a} such that 
{a} U {a} = {u}. There are 2™i ways to partitionate {u} into {a} and {a}, and each possibihty 
is simply recasted as 

{a} = {ui, . . . ,u,ii} and {a} = {ui, . . . , J, (54) 
with mi = ni + ni. For each partitioning we also define the functions 



- /_ \-\Li + l ni ni I - \\ 



and following |15j the structure constant for three-point functions at leading order in planar A/" = 4 
SYM is given by 



under the constraint mi = m2 + va-^. The normalisation factors are 

mi mi mi 



n : : : ^> i(K},M), 

m2 r7l2 ^'^2 



xii n : "i^^V l.(k},w), 

1 ■ ii ~ + ^) 

ma ma ma 

AAj =(-1)™^ n H^,r'''aiw*)-^^aiw,)b{w*) J] H 



a(u'i — Wj + ^) 



j=i i=ij=z+i "-y^^^ ^3 ' 2) 

xR n T^^^^a W) , (57) 

where A^^ denotes the scalar product of Bethe vectors (|46| for a lattice of length L after the limit 
(|48l) . The function A''^ turns out to be explicitly given by 



1 " 

^(M,{A))=— 



(27ri)2" ^11^ a(Ai - i^j + ^)a{i^j - X^ + ^)a{u, - i/j + ^)2 

n?"i n?=j+i - + §)6(wj - Wi + ^) 



diDi . . . dw 



2n- 



u-Zi n;=i bi^^ A, + i)6K - 1., + i) 



n n n 

X JJ[a(wi)fo(wn+j)]^ a(Wi - + ^)a(^'j - Wn+i + ^) 

2—1 i—1 j—1 

n n 

X ]^ b{iyj - Wi + ^)b{w2n+i-i ~ + \)a(v^ - Wj + ^)a{w2n+i-j - Vi + ^), 

i=i j-i+i 

(58) 

with the functions a and b as described in (|50p . The number of excitations n is understood by the 
number of rapidities in the argument. 
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In its turn the function 6*2(0!) is given by 



"1 "12 "12 "12 I / _ J_ j_\ 

(59) 



while 



^3(a)-n^K)"'"'«(^^)n«(^^-)"'"-^^(^^-)n n ) \ \t 

j=i 3=1 '^3 + 2) 

X ^ni3({^^})-^ni3({^}) ^ n 1,7713 ■ 

(60) 

In the expressions (|59p and (|60p . (5^.^ stands for the standard Kronecker delta and consists 
of the six-vertex model partition function with domain wall boundaries (|36p for a lattice with 
dimensions 7713 x 7713 in the limit (|48 p -(|50 p previously discussed. For completeness it is given by 



"13 

X n Hw^y-^a{w,^^-' dwi... dwm, , (61) 



recalling that a and b correspond to the functions in ((50)) . 

In this way, to compute the structure constant leading term 0^23 according to the prescriptions 
of [TS], one needs to substitute (|58l) in (|57)) and ([55]) . and then insert the results into ([55)1 . Similarly 
we also need to substitute (|5T|) in (150]) and replace the results in (|56p . In summary, the evaluation 
of CJ23 is th^s reduced to the evaluation of the countour integrals and (|5T|) . 



7 Conclusion 

The main result of this paper is to write the partition function for a general Z-invariant six- 
vertex model as a multiple contour integral over a factorised polynomial kernel. This allows us 
to write a new representation for the partition function of the six vertex model with domain 
wall boundary conditions as a multiple contour integral. Likewise we derive an off-shell multiple 
integral expression for the scalar product of two Bethe vectors of the six-vertex model transfer 
matrix. 

The study of two-point functions mAf — A Super Yang-Mills (SYM) has advanced dramatically 
due to the presence of integrable structures, and the computation of three-point functions at weak 
coupling seems to be benefited as well [TS]. In this sense we have also illustrated in Section IH] how 
the integral formulaepresented in Sections 15.11 and 15.21 can be used in that context. Although in a 
different fashion, recently the connection was realised between the partition function of a system 
of statistical mechanics and the computation of norms of Bethe vectors in a particular limit |30] . 
We hope the results presented here to help understanding this relation in more general cases. 

The integral formulaewe obtain are very similar to those in multi-variate polynomial solutions of 
the q-deformed Knizhnik-Zamolodchikov equation, a connection which deserves future exploration. 
In |19fl23j methods have been developed to compute with such multi-variate expressions which we 
hope will be applicable in the current context as well. We further hope that our representation 
for the scalar product can be extended to compute correlation functions, providing an alternative 
to the methods developed in , and also constitutes an approach which may be generalised 

to higher rank solvable lattice models. 
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A Details of scalar product computation 

We consider the result of Baxter [Tl], specialised to a 2n by L square lattice, i.e. we have 2n + L 
lines and 4n + 2L rapidities, (see Fig. [6] ) 

Thus we have the set of right rapidities (in a clockwise direction), using the notation = 

^J■L-^+l 

U = {ui,...,w2„+l} = {/ii,...,^L,A„,...,Ai,z/i,...,r/n} (62) 
and the full set of rapidities 

V = {Ui, . . . ,M2n+L} 

= {^i,...,/XL,A„,...,Ai,i/i,...,r/„,/2i -7, -7,j^i -7,---,i^i -7,Ai -7, ...,A„ -7} 

(63) 

We have the out arrows appearing in three distinct blocks. The positions of the out arrows are 

X = {1, . . . , L, L + n + 1, . . . , L + 2n,2L + 2n + I, . . . ,2L + 3n} (64) 

With these specialisations, we can see that a number of terms of the original {L+2n)\ permutations 
will vanish. 

If we first consider the (t>{wi, Xi) terms, we find that the first L of these (corresponding to the 
first set of out arrows, associated to the /i rapidities), will be zero, except for the term arising 
from the poles at Wi = fXi for i = 1, . . . , L. It is then easy to perform the first L integrals, and we 
find that the integral formula for the scalar product becomes: 

(_l)(i-i)»(sinh7)2" A 



(27ri)2n ^11^ a(A, - Vj)a{vj - Xi)a{i^i - Vjf 

j' ■■■ j' dl«L+l • • • dWL+2„ ^ \\ F{wL+i,WL+n+i)^ 



rii"! n?=»+l a{wL+j - WL+i)b{wL+j - WL+i) 
X 

n^=i lVj=i b{wL+t - Xj)b{wL+i - lyj) 

n L n n 

X ]^ a{wL+i - ^J.j)b(p■j - WL+n+i) Y\_ II a(wL+-t - l'j)a{vj - WL+n+i) 
i=l j—1 i—1 j—1 

n n 

X II II b(l^j - WL+i)b{wL+2n+l-i - '^j)a{iyi - ■WL+j)a{wL+2n+l~j - I'i) ■ 

(65) 



with 



^ TT a{wL+i - fij)a{wL+n+i - l^j)a{P] - WL+i)a{^ijWL+n+i) 

^ TT- aj-WL+i - Xj)a{wL+n+i - Aj)a(Aj - WL+i)aiXj - WL+n+t) ^gg^ 
jj[ (a(A» - Aj))2a(Aj - l'J)a{|y^ - Xj) 
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If we now consider the 4'{wi,Xi) terms corresponding to the next two sets of out arrows, we 
find that we have terms that cancel all poles at Vj except for j = i for WL+i and that also the only 
poles associated to the WL+n+i variables remaining in the set of v rapidities are again those when 
i = j. 

In summary, we have poles at: 

^ Hi for i = 1, . . . , L, 

wl+i ^ i^i or X^. for z = (67) 

WL+n+i = i^n+i-i or A^„^i_i for i = 1, . . . , n. 

This enables us to further simplify (|65p . It is easy to see that the factor F is equal to one for every 
non- vanishing choice of poles in ((67)) . and thus we obtain the expression (|46l) . 
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